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We consider the spin-1/2 antiferromagnetic Heisenberg chain in a staggered magnetic field describ- 
ing materials such as copper benzoate and copper pyrimidine dinitrate. Using the density-matrix 
renormalization group for transfer matrices (TMRG) we calculate the magnetization of these materi- 
als at finite temperature and arbitrary magnetic field. These results are in excellent agreement with 
experimental data, allowing for a determination of the inhomogeneity parameter c of copper ben- 
zoate (c = 0.043) and copper pyrimidine dinitrate (c = 0.11). The TMRG approach can be applied 
to rather low temperatures yielding singular field and temperature dependences of susceptibilities. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.30.Gw, 75.50.Xx 



One-dimensional quantum magnets have been of theo- 
retical and experimental interest in recent years because 
of the rich variety of different magnetic ground states, 
such as quantum critical behavior or gaps in the spin exci- 
tation spectraii^i 3 .. The ideal spin-1/2 antiferromagnetic 
Heisenberg chain (5=1/2 AFHC) with uniform nearest 
neighbor exchange coupling is of particular interest, since 
it is exactly solvable using the Bethe ansatz (BA)^*^. Its 
ground state is a spin singlet with gapless excitations and 
quasi-long ranged ground state correlations, hence even 
small perturbations can change the physical properties 
fundamentally. 

In experimental realizations of spin chains like copper 
pyrimidine dinitrate [PM Cu(N0 3 )2(H 2 0)2] n (CuPM) 
and copper benzoate Cu(C6H 5 COO)2 • 3H 2 additional 
terms in the Hamiltonian result from the lack of inver- 
sion symmetry. As a consequence of the residual spin- 
orbit coupling the Dzyaloshinskii-Moriya (DM) interac- 
tion and an alternating g tensor have to be taken into 
account^. This gives rise to an effective staggered field 
h s perpendicular to the applied magnetic field H. The 
Hamiltonian is written &s£& 
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with the coupling constant J, the effective uniform field 
h u = gfisH/J and the induced effective staggered field 
h s = ch u . For a given field axis the effective g and the 
inhomogeneity parameter c are determined from the DM 
interaction D and the alternating g tensor g = g u ± g s 
in the following way^ g = g u H / H and 
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In essence, for materials like CuPM and Cu-bcnzoate, 
the presence of DM interaction and staggered g tensor 
induces new symmetry breakings with respect to the 



magnetic properties. In terms of principal axes, conven- 
tionally it is distinguished between the crystallographic 
unit cell with the axes abc, the coordinate frame a'bc', 
where the uniform part of the g tensor is diagonal, and 
the system of the principal magnetic axes a"bc"£&. For 
both materials the g tensor and the DM vector have 
been derived previously, as have the different coordinate 
frames^. Hence, these quantum spin systems are per- 
fectly suited as model compounds for a quantitatively ex- 
act comparison between theory and experiment. Corre- 
spondingly, in recent years various experimental and the- 
oretical studies on static or dynamic properties of stag- 
gered S=l/2 AFHC have been carried ou1*i2*ii*i£, even 
though data analysis was essentially limited to zero tem- 
perature. Therefore, in this combined theoretical and 
experimental study we provide a detailed investigation 
of the effect of temperature on thermodynamic proper- 
ties of staggered S=l/2 AFHCs. For the first time, we 
derive the induced staggered magnetization for the stag- 
gered S=l/2 AFHC from the low temperature regime 
ksT <C J up to the paramagnetic range ksT > J in fields 
up to saturation. Since the effective g factor and the c 
factor depend on the orientation of the magnetic field H 
with respect to the crystal axis @, we have to calculate 
these physical parameters for each experimental config- 
uration separately. For this reason, here we calculate 
the magnetization for different crystal alignments by use 
of the transfer matrix renormalization group (TMRG) 
method and compare these data with experimental data 
of CuPMi 3 . and Cu-benzoate. The good agreement of 
numerical and experimental data allows for the determi- 
nation of c and g factors. 

To study thermodynamical properties at finite tem- 
perature the TMRG method provides a powerful numer- 
ical tool, because the thermodynamic limit is performed 
exactly and in contrast to quantum Monte Carlo tech- 
niques there is no "minus sign" problem. The idea of the 
TMRG method is to express the partition function Z of 
a one-dimensional quantum system by that of an equiva- 
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lent two-dimensional classical model by a Trotter-Suzuki 
mappingi^^. For the classical model a suitable quan- 
tum transfer-matrix (QTM) can be defined which allows 
for the calculation of all thermodynamical quantities in 
the thermodynamic limit by considering solely the largest 
eigenvalue of this QTM. Here we use the Trotter-Suzuki 

mapping^ yielding Z = lim M ^oo Tr { [Ti(e)T 2 (e)] M/2 } , 

where T li2 (e) = T RX exp [-eH + 0{e 2 )] , e = fi/M with 
(3 being the inverse temperature and M a large inte- 
ger number. Tr l denote the right- and leftshift oper- 
ators, respectively. The bulk free energy is given by the 
largest eigenvalue Ao, / = — Tin Ao . Expectation values 
of local operators, like the homogeneous and the stag- 
gered magnetization, can be expressed in terms of the 
left and right eigenvectors belonging to the largest eigen- 
value. For achieving lower temperatures the length of the 
transfer matrix is increased in imaginary time direction 
by an application of the infinite density matrix renor- 
malization group (DMRG) algorithmic. In all following 
calculations we have retained between 64 and 100 states 
and used e = 0.05. There are no finite-size effects as 
the system size is strictly infinite, however numerical er- 
rors due to the discretization of the imaginary time axis 
and the DMRG truncation are generally within the line 
width used in our plots as has been checked by compar- 
ison with the exact BA data (sec Fig. For details of 
the algorithm the reader is referred to Refs>i§iiLi2ii2i2ii. 

In the following we compare our numerical results 
with high-field magnetization data measured at the Lab- 
oratoire National des Champs Magnetiques Pulses in 
Toulouse in pulsed magnetic fields up to /io-ff = 53 T 
(CuPM) and 38 T (Cu-benzoate), respectively (experi- 
mental details in Refii^), and with results obtained by 
means of 13 C NMR2I. In the presence of a staggered 
magnetic field h s = c h u induced by a uniform field 
h u , the magnetization m is given by the superposition 
of the uniform m u and the staggered m s magnetization 
component 
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m u +cm s 



(3) 



where c is the above inhomogeneity parameter. The mag- 
netization m can be measured experimentally. On theory 
side we can determine the magnetization m by using © 
and calculating the uniform m u and the staggered mag- 
netization m s separately by TMRG methods. 

Recently, CuPM has been identified as a 5=1/2 
AFHC, with a magnetic exchange parameter J/ks = 
(36 ± 0.5)KS (in the following we use J/fcs = 36. 5K). 
In the uniform 5=1/2 AFHC model the (magnetic) sat- 
uration field is calculated according to the formula^ 



H c = AJS/gfiB 



(4) 



The g tensor and the DM vector of CuPM have been 
derived from ESR measurements and single-crystal sus- 
ceptibility studies^. In the coordinate frame a'bc' the g 
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FIG. 1: Plot of theoretical (solid line) and experimental (open 
circles) magnetization data for CuPM at T = 1.6 K along the 
a" direction. Dashed line shows TMRG data with a misalign- 
ment of 5° at T = 1.6 K. Note that deviations are strongest 
in the vicinity of the saturation field. Lower plot: difference 
between TMRG and BA data. 



tensor takes the form 



9 = 





2.149 ±0.127 
±0.127 2.287 



= 9u ± 9s 



and the DM interaction vector reads 



D = 0.139J(-0.4115, 0,0.9114) 



(5) 
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If the magnetic field H is directed along the a" direc- 
tion of the CuPM crystal, it will behave like an ideal 
5=1/2 AFHC (c = 0). In this case, from © we obtain 
the effective g factor as g = 2.14, with saturation field 
HqH c = 50.7 T Q. In Fig. 2] we present the magne- 
tization curve m of CuPM as function of the magnetic 
field at T = 1.6 K along the a" axis, as well as the differ- 
ence between TMRG and the exact BA results^. Since 
this difference is always smaller than 2 • 10 -3 the TMRG 
results describe the exact BA results very well. Com- 
paring the experimental data with the TMRG results we 
find good agreement with the uniform 5=1/2 AFHC for 
magnetic fields up to /.iqH < 45 T. For higher fields the 
Dzyaloshinskii-Moriya interaction and the alternating g 
tensor may take effect on small crystal misalignments, 
so that CuPM does not behave like a uniform 5=1/2 
AFHC. Assuming a misalignment of 5° of the magnetic 
field with respect to the a" axisi 3 . the physical param- 
eters have to be changed to J/ks = 36.5 K, g = 2.13 
and c = 0.01 J5J). With these slight modifications the 
calculated magnetization curve of CuPM describes the 
experimental data well, see Fig. ^ 

Along the c" axis the effect due to the induced stag- 
gered field is largest. According to& a gap is induced 
A cx /i 2 / 3 with multiplicative logarithmic corrections. At 
low T the magnetization receives strong contributions 
from the staggered component with singular behavior. 
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FIG. 2: Susceptibility \ °f CuPM as a function of H at 
T = 1.6 K calculated by TMRG. Solid (dashed) line: suscep- 
tibility with magnetic field along c (a") axis. Inset: x °f 
CuPM as a function of T at ft = calculated by TMRG. 
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FIG. 3: Theoretical and experimental magnetization curves 
for CuPM at T = 1.6 K along the c" direction. Solid line: 
magnetization m calculated for c = 0.11, dashed line: uni- 
form magnetization m u , dotted line: staggered magnetiza- 
tion m s , open circles: experimental data. Inset (a) depicts 
the magnetization for small fields calculated by TMRG and 
ED (dashed-dotted line) in comparison to the experimental 
data, inset (b) displays the temperature dependent staggered 
magnetization of CuPM for a magnetic field fj,oH = 9.3 T 
applied along the chain axis (solid line: TMRG results; open 
circles: experimental data). 



For T = the dependence m s oc h 1 / 3 was found and for 
h = at low T the susceptibility x(T) oc 1/T with mul- 
tiplicative logarithmic corrections. We stress that this 
behavior is the result of strong correlations despite vague 
similarities with paramagnetic impurities. Susceptibility 
data are shown in Fig. [21 Note however, the direct appli- 
cation of DMRGi£ to the Hciscnberg chain with full DM- 
terms does not give evidences of logarithmic corrections. 
Along this axis we obtain the parameters c = 0.11 and 
g = 2.19 with saturation field /io-ffc = 49.6 T according to 
. Figure [3]shows magnetization curves for CuPM with 
H along the c" axis at temperature T = 1.6 K. For all 



magnetic fields the magnetization calculated by TMRG is 
in good agreement with the experimental data. In partic- 
ular at small fields, in comparison with the results of ex- 
act diagonalization (ED)A2i the TMRG describes the ex- 
periment much more accurately, because it has no finite- 
size effects by construction (inset (a) of Fig. |2J. We 
stress that the ED results were calculated for T = 0, 
whereas the TMRG results are calculated at the same 
finite temperatures chosen in the experiment. 

For magnetic field H parallel to the chain axisS of a 
CuPM crystal we find the parameters c = 0.083 and 
g = 2.117. In inset (b) of Fig. [3] we present the ex- 
cellent agreement of the calculated and the measured^! 
staggered magnetization curves for magnetic field fj,oH = 
9.3 T. Since these results are in good agreement with the 
experimental data, wc verified the physical parameters 
for CuPM, especially c = 0.11 along the c" axis with 
accurancy ±0.01. Note that this is consistent with re- 
sults based on EDi£, but rules out c=0.08 obtained in 
[23] from ESR. 

Cu-benzoate is another example for a staggered S = l/2 
AFHC, with a behavior very similar to CuPM. Only, for 
Cu-benzoate the coupling constant is smaller by a fac- 
tor - 2 than for CuPM, that is J/k B = (19 ± 0.5) Ki£. 
In consequence, for this material, in a temperature de- 
pendent high magnetic field study, the regimes from the 
fully staggered one fcgT <C J up to the paramagnetic 
range fcgT > J in fields up to saturation are accessible 
and can be compared to the results from TMRG calcula- 
tions. From ESR measurements the g tensor in the a"bc" 
coordinate system takes the form: 

2.115 ±0.0190 0.0906 
±0.0190 2.059 ±0.0495 ) (7) 
0.0906 ±0.0495 2.316 

Moreover the DM interaction has been determined^ 

D = J (0.13, 0, 0.02) . (8) 

Cu-benzoate will behave like a homogeneous 5=1/2 
AFHC (c = 0) with a coupling constant J/fcs = 19 K, 
if the magnetic field is along the a" axis. The saturation 
field is calculated to \i§H c = 26.5 T by using g = 2.13 
corresponding to this axis. In Fig. 0]we show the magne- 
tization ® as a function of the magnetic field at differ- 
ent temperatures. Overall, the TMRG results agree very 
well with the experimental data obtained in our high- 
field magnetization study, especially with respect to the 
saturation magnetization. The field induced gap of Cu- 
benzoate is largest for magnetic field along the c" axis. In 
this case we use J/kB = 18.9 K and calculate the effective 
g factor to g = 2.32, with a corresponding saturation field 
HqH c = 24.2 T. By using (J2J) the c factor is determined to 
c = 0.043. The inset of Fig. 0] shows the magnetization 
m (J3J) of Cu-benzoate along the c" axis at different tem- 
peratures. Altogether, within experimental resolution we 
find that the magnetization of Cu-benzoate is nicely de- 
scribed by TMRG results using c = 0.043±0.01 (covering 



4 




FIG. 4: Homogeneous magnetization of Cu-benzoate along 
a" axis at different temperatures (1.7K, 4.2K, 12K and 25K). 
Solid (dashed) lines show TMRG (experimental) data. Inset: 
Magnetization of Cu-benzoate along c" axis at different tem- 
peratures (1.7 K, 4.2 K, 12 K and 25 K). The dotted lines 
show the magnetization at T = 12. 8K and T = 22. 5K calcu- 
lated by TMRG. 

the value c = 0.034 resulting from the DM vector used in 
RefiiS). The somewhat larger deviations between exper- 
iment and theory at higher temperatures (T = 25 K) can 
be ascribed to a larger experimental uncertainty (±3 K) 
in this temperature range. 

In Fig. [3] it can be clearly seen that in the low tem- 
perature staggered regime (1.7 K = 0.09fcsT/J) there 
is a prominent staggered component, with a maximum 
value of about 0.6 /ib/Cu atom at 20 T. This staggered 
magnetization is gradually wiped out as temperature is 



increased up to 25 K = l.3ksT/ J, viz., the conventional 
paramagnetic regime. 

In summary, a comparative analysis of numerical 
TMRG and experimental data for Cu-benzoate and 
CuPM at low temperatures was presented. In our com- 
parison we found very good agreement between our nu- 
merical and experimental data. The field dependence of 
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FIG. 5: Staggered magnetization of Cu-benzoate along c" axis 
at different temperatures (1.7 K, 4.2 K, 12 K and 25 K). 

the magnetization was found to show interesting charac- 
teristics at high fields corresponding to saturation and fol- 
low field due to quantum criticality of the 5=1/2 AFHC. 
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